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1.0  INTRODUCTION 


The  buckling  analysis  of  shells  has  in  recent  years  gained 
considerable  interest.  The  current  demand  for  lightweight 
efficient  structures  requires  the  engineer  to  closely  examine  the 
response  of  the  shell,  particularly  for  those  states  of  stress 
that  lead  to  buckling.  It  has  been  observed  that  as  the  design 
of  a  shell  structure  is  optimized,  the  structure  becomes  more 
prone  to  instability.  Further,  small  perturbations  in  load  or 
geometry  will  not  permit  the  structure  to  sustain  its  load,  and 
failure  will  occur.  In  some  cases,  linear  methods  did  not  give 
satisfactory  results  and  nonlinear  theories  were  employed  in  an 
attempt  to  improve  the  correlation  between  prediction  and 
experiment.  This  led  to  postbuckling  and  imperfection  studies, 
which  in  turn  provided  a  measure  of  the  load-carrying  capacity  of 
the  structure.  Unfortunately,  these  analyses  did  not  give  an 
accurate  prediction  of  the  residual  strength  of  the  shell,  but 
did  give  an  adequate  assessment  of  its  behavior  just  prior  to 
buckling. 

The  thin  stiffened  and  unstiffened  circular  cylindrical 
shell  has  been  extensively  used  in  the  aerospace  and  energy 
industries.  The  cylinder  was  found  to  be  particularly 
imperfection-sensitive  to  axial  compression.  This  sensitivity 
was  attributed  to  the  deviation  of  the  shell's  geometry  from  a 
"perfect"  shell.  This  deviation  from  the  true  or  perfect  shape 
is  termed  an  imperfection.  Thus,  the  phrase  "imperfection 
sensitivity"  has  come  into  popular  use. 

During  the  initial  phase  of  this  research,  a  review  of  the 
literature  was  conducted  and  a  report  prepared  [l].  In  this 
review,  the  analysis  methods  derived  by  Donnell  and  Koiter  were 
closely  examined.  Donnell  derived  a  set  of  nonlinear  equations 
that  govern  the  response  of  an  axially  loaded  cylinder.  Koiter 
examined  the  behavior  of  the  shell  at  adjacent  states  of  stress 
to  the  buckling  load  and  hypothesized  the  buckling  character¬ 
istics  through  these  states.  Koiter's  theory  is  based  on  the 
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initial  postbuckling  behavior  of  a  perfect  structure.  The  theory 
leads  to  an  equation  for  a  load  versus  postbuckling  displacements 

V  A*  =  1  +  a  e+  b  ^  j  ^ 

in  which  x/ X*  is  the  ratio  of  applied  load  to  the  classical 
buckling  load,  e  is  the  amplitude  of  displacement  of  the  buckled 
mode,  and  a  and  b  are  coefficients  that  depend  on  the  geometry  of 
the  structure,  prebuckling  behavior,  and  type  of  loading.  Equa¬ 
tion  (1)  applies  to  cases  with  unique  eigenvalues. 

The  coefficient  b  is  more  Important  than  the  coefficient  a 
in  most  cases  Involving  thin  shells;  therefore,  the  coefficient  a 
usually  takes  on  a  value  of  zero,  particularly  for  shells  of 
revolution.  Often  referred  to  as  Koiter's  parameter,  the 
coefficient  b  essentially  determines  the  Imperfection  sensitivity 
characteristics  of  the  shell.  When  b  is  less  than  zero,  X/X*  is 
less  than  unity,  thus  predicting  a  reduction  in  load-carrying 
capacity  from  that  of  the  perfect  shell.  Computer  methods  are 
often  required  for  the  evaluation  of  the  imperfection  sensitivity 
parameter  b.  Only  one  computer  program  system,  written  by  Cohen 
[2],  is  specifically  designed  to  examine  the  imperfection  sensi¬ 
tivity  of  general  shells  of  revolution.  In  using  this  program, 
one  first  has  to  execute  a  static  program  module  to  determine  the 
prestress  state.  These  results  are  passed  to  a  stability  program 
which  uses  a  subroutine  to  establish  the  imperfection  sensitivity 
of  the  shell.  Although  Cohen's  procedure  appears  to  be  correct, 
it  has  not  been  widely  accepted  by  the  engineering  community 
because  the  program  is  not  user-friendly.  There  remains  a  need 
for  an  efficient  computer  procedure  that  is  easy  to  use  and  yet 
suitable  for  preliminary  design. 

In  the  following  sections  the  theoretical  basis  of  four  such 
computer  programs  are  presented.  The  programs  are  designed  to 
evaluate  the  buckling  characteristics  of  a  cylindrical  shell 
stiffened  by  rings  or  by  rings  and  stringers.  Equivalent  ortho¬ 
tropic  membrane  and  bending  stiffnesses  are  derived  through  a 


model  in  which  the  ring  and  stringer  stiffeners  are  "smeared" 
over  the  shell.  For  light  and  medium  stiffening,  where  the 
cross-section  area  of  the  stiffener  divided  by  the  stiffener 
spacing  is  less  than  one-half  the  thickness  of  the  shell,  the 
smeared  model  has  been  found  to  be  adequate. 

For  ease  of  identification,  each  computer  procedure  has  the 
project  name  "PVRC"  followed  by  the  initial  of  the  principal 
contributor : 


Boros  [3] 
Hutchinson  [4] 
Khot  [5] 

Arbocz  [6] 


PVRCB 

PVRCH 

PVRCK 

PVRCA 


The  capabilities  of  these  programs  are  listed  in  Table  1 . 

Because  of  the  applicable  theories  employed,  only  PVRCH  considers 
the  boundary  conditions  at  the  ends  of  the  cylinder.  For  the 
remaining  procedures,  simply-supported  boundary  conditions  are 
approximated . 

All  four  programs  are  based  on  Donnell's  equations,  but  dif¬ 
ferent  solution  procedures  are  employed,  and  different  loading, 
imperfection  geometry,  and  boundary  conditions  allowed.  Each  of 
these  procedures  will  be  discussed  in  detail. 

The  first  program,  PVRCB,  examines  the  buckling  character¬ 
istics  of  a  constant  thickness  circular  cylinder  stiffened  by 
rings  only.  The  imperfection  is  assumed  to  be  axisymmetric  and 
the  effects  of  the  boundary  condition  are  not  considered.  A 
secondary  option  is  included  for  a  two-mode  imperfection  (one 
axisymmetric  and  the  other  nonaxisymmetric) .  The  applied  load 
can  only  be  an  axial  stress. 

The  second  program,  PVRCH,  examines  the  buckling  strength  of 
a  circular  cylinder  stiffened  by  stringers  and  rings.  The  shell 
can  have  either  clamped  or  simply  supported  conditions.  The 
loading  can  be  either  an  axial  stress  or  normal  pressure.  Only 
an  axisymmetric  imperfection  geometry  can  be  considered. 

The  third  program,  PVRCK,  also  considers  a  ring  and  stringer 
stiffened  circular  cylinder;  however,  the  cylinder  is  comprised 
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TABLE  1.  IMPERFECTION  SENSITIVITY  PROGRAMS 


Capability  Parameters 

Program 

Rings 

Stringers 

;  Boundary 
Conditions 

i  Layer 
(no./type) 

Loading 

Imperfection 

Geometry 

PVRCB 

X 

1/ISO 

AX 

symmetric 

nonsymmetric 

PVRCH 

X 

X 

SS,  CL 

1/ISO 

AX ,  P 

symmetric 

PVRCK 

X 

X 

SS 

5 /ORTHO 

AX , P  ,T 

nonsymmetric 

PVRCA 

X 

X 

~  SS 

1/ISO 

AX 

arbitrary 

Legend 


AX 

axial  load 

P 

pressure 

T 

torsion 

ISO 

isotropic 

ORTHO 

orthotropic 

SS 

simply  supported 

CL 

clamped 

4 


of  a  layered  orthotropic  material.  Boundary  conditions  are  not 
imposed.  The  loading  can  be  any  combination  of  axial  compres¬ 
sion,  pressure  and  torsion.  An  imperfection  amplitude  acting  in 
each  critical  longitudinal  and  circumferential  response  mode  is 
assumed  and  the  buckling  load  is  computed. 

The  fourth  program,  PVRCA,  uses  Koiter's  method  for 
calculating  the  buckling  load  of  axially  compressed  stiffened 
cylindrical  shells  with  a  given  asymmetric  imperfection. 

Boundary  conditions  are  ignored.  No  other  loads  are  considered. 

Fortran  listings  of  the  four  computer  programs  are  included 
as  appendixes  (see  Volume  II  of  this  report).  In  addition, 
common  driver  subroutines  are  listed  in  Appendix  E.  These 
programs  are  designed  to  execute  in  an  interactive  mode. 
Therefore,  no  preparation  of  an  input  data  file  is  required.  The 
user  is  automatically  prompted.  Copies  of  the  CRT  displays  are 
provided  to  illustrate  the  use  of  each  program. 

The  governing  equations  of  a  cylindrical  shell  of  revolution 
are  presented  in  Section  2.1.  The  four  sections  that  subse¬ 
quently  follow  are  the  detailed  derivations  that  form  the  basis 
of  the  solution  procedures  for  PVRCB,  PVRCH,  PVRCK,  and  PVRCA, 
respectively.  The  significance  of  certain  assumptions  as  to 
loading,  boundary  conditions,  construction  of  the  shell's  wall, 
etc.  will  also  be  discussed  at  the  appropriate  point  in  the 
development . 

Numerical  examples  obtained  from  the  first  three  developed 
computer  programs  for  four  ring  stiffened  and  four  stringer 
stiffened  configurations  are  presented  in  Section  3.0.  A 
separate  verification  example  is  also  presented  for  the  program 
PVRCA.  Section  4.0  discusses  the  behavior  of  the  Imperfection 
sensitivity  parameter  and  the  applicability  and  usefulness  of  the 
four  programs.  Finally,  conclusions  and  some  observations  are 
presented . 
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2.0  TECHNICAL  DISCUSSION 
2 . 1  Governing  Equations 

Consider  a  cylindrical  shell  of  radius  R,  length  L,  and 
thickness,  t.  If  the  usual  assumption  of  thinness  is  made 
((t/R)  <<  1),  then  the  governing  equations  of  thin  shell  theory 

can  he  developed.  For  the  coordinate  system  shown  in  Figure  1, 
approximate  nonlinear  strain-displacement  and  curvature- 
displacement  equations  can  he  written  as: 


i  e  }  - 


1  U,x+?<\*> 


<  e 


2 


xy 


=  V  +  i.  (w  )2  _  X 
,y  2  U,y'  R 


U  +  V  +  w  w 


{<}  = 


¥ 


=  W 


,  xx 


K 

xy. 


»yy 
2  w 


»xy 


(2) 


in  which  U,  V,  W  are  the  axial,  circumferential  and  radial 
displacements  of  the  reference  surface,  and  the  comma  denotes 
differentiation  with  respect  to  the  subscripted  variable. 

The  equilibrium  equations  derived  by  Donnell  are: 


N  +  N  =0 

x, x  xy,y 

N  +  N  =0 

y , y  xy , x 


•s-N+NW  +2N  W  +NW 
R  y  ™ x  , xx  +  xy  w,xy  +  y  ,yy 


(3) 


+  M  +  2  M  +  M  =0 

x , xx  xy ,xy  y,yy 


where  N  =  {N  ,N  ,N  (the  stress  resultants)  and  M  = 

r  i  TX  /  ^  \ 

|M  ,M  ,M  }  (the  moment  resultants).  These  resultants  are 
1  x ’  y ’  xy  J 

usually  defined  along  the  reference  surface  of  the  shell  and  are 
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Reference  Surface 


Reference 


Figure  1 


(b) 


Geometry  of  the  Shell 
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the  average  and  linear  variation  of  the  inplane  stresses  through 
the  thickness  of  the  shell,  respectively. 

The  constitutive  equations  for  a  layered  orthotropic  shell 
of  revolution  can  be  written  in  matrix  form  as: 


where  the  membrane,  coupling,  and  bending  stiffness  of  the 
composite  shell  can  be  defined,  respectively,  as: 


N 


(k) 


=  Ji  V  '  (h 


k+1  “  hk^ 


N 


K 


13  ■  7  X  Bi.i(k)  (hk+l  -  hk> 


k=1  xJ 


N 


D 


00 


13 -y  X  Bi3W(hk+i-^ 

(k) 

B..v  =  compliance  moduli  in  each  layer  k 


id 


(e.g.,  B 


E 


1 1 


1  1  1  -v1  v2 


) 


h^  =  distance  from  reference  surface  to  each  layer  k 


The  matrix  coefficients  C. .,  K. .,  and  D. .  can  be  modified  to 

^■3  i  0 

include  the  effects  of  the  stiffening  element: 


c.  . 

1J 

Cij 

+ 

6.  . 
1J 

EA 

b 

K.  . 
1J 

Kij 

+ 

£ 

EAe 

0  .  . 
10 

b 

/V 

DiJ 

Dij 

6  .  . 
1J 

El 

(5) 
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where 


=  Kronecker  delta  function 
E  =  modulus  of  elasticity  of  the  stiffener 
A  =  cross-sectional  area  of  the  stiffener 
b  =  distance  between  stiffeners 

e  =  eccentricity  of  the  stiffener 

I  =  moment  of  inertia  of  the  stiffener 

If  i=g'  =  1  ,  the  above  properties  are  for  stringers;  if  i=j=2,  they 
represent  ring  properties  and  both  are  "smeared"  over  the  span  b. 
When  i*j  the  Poisson's  effect  of  a  shell  should  be  considered. 
Plugge  [7]  develops  an  anisotropic  approximation  for  a  stiffened 
shell  but  does  not  include  Poisson's  effect.  Hutchinson  and 
Araazigo  [8],  however,  developed  equations  which  account  for  a 
partial  effect. 

In  order  to  examine  the  imperfection  sensitivity  of  a  shell, 
the  imperfection  geometry  must  be  introduced  into  Eqs.  (2)  and 
(3).  The  radial  variation  from  a  perfect  form  is  usually  taken 
as  the  imperfection,  W.  The  total  radial  displacement  in 
Eqs.  (2)  and  (3)  is  therefore  comprised  of  a  response  and  an 
imperfection.  Thus,  additional  terms  must  be  included  in  the 
above  relationships  so  that  there  is  no  stress  when  at  zero  load. 
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2.2  Program  PVRCB 

This  program  examines  the  buckling  characteristics  of  a 
constant  thickness  isotropic  circular  cylinder.  Only  ring 
stiffeners  are  considered,  and  their  properties  are  smeared  as 
described  in  the  previous  section.  The  out-of-plane  bending 
stiffness  and  the  torsional  rigidity  of  the  ring  are  neglected. 
The  applied  load  can  only  be  an  axial  stress  (a),  although  an 
applied  pressure  (p)  will  be  carried  during  the  derivation. 
Boundary  conditions  are  not  explicitly  stated,  but  simply 
supported  conditions  are  implied. 

The  development  of  the  governing  equations  and  the  basic 
solution  procedure  is  based  upon  the  works  of  Tennyson,  his  stu¬ 
dents  [3],  and  colleagues  [9]«  Boros  [3]  presents  the  derivation 
for  the  effects  of  the  axisymmetric  imperfection,  and  Hutchinson 
[9]  presents  the  effects  of  a  nonaxisymmetric  imperfection. 

Assume  an  Airy's  stress  function  F,  defined  as: 


N 


x 


F 

>yy 


N 

y 


F 

,xx 


N  =  -F 
xy  ,xy 


(6) 


The  Donnell-von  Karman  shell  equations  [10]  take  the  form,  for 
compatibility: 


L„[F]  -  Ln[w]  -  ¥  +  ¥  ¥  +  ¥  ¥  +W 

HL  Q  ,xy  ,xx  ,yy  o,xx  ,yy 


Q 

A 

-  2  ¥ 


¥  n 
o ,  xy  ,xy  =  0 


¥ 

0 ,yy  ,xx 
(7) 


and  for  equilibrium: 


L  [¥]  +  L_[F]  -F  ¥  -F  ¥  +2F  ¥ 

D  Q  ,xx  ,yy  , yy  ,xx  , xy  , xy 

A  A  A 

-  F>xx  ¥0)yy  -  F,yy  ¥0,XX  +  2  F>Xy  W0,xy  +  p  =  0  (8) 
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where 


The  formulas  for  effective  stretching  and  bending  stiffnesses  are 
given  in  Ref.  [8]  and  repeated  herein  in  Table  2. 

2.2.1  Axisymmetric  Imperfection 

Assume  an  axisymmetric  imperfection  of  the  form: 


(10) 


S  t  cos 


o 


where 


k  =  1 ,2,3 


L  =  length  of  the  shell 
=  imperfection  amplitude 
t  =  thickness  of  the  shell 

Prior  to  buckling,  the  radial  deflection  and  stress  function  are 
assumed  to  be  of  the  form: 


F 


W 


^  oR  +  W* 

1  2 

-  oty  +  F* 


(11) 


where 


a 


axial  stress 


W*  =  B  cos 


2p  x 
o 


(12) 
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TABLE  2 


DEFINITION  OF  STRESS  COEFFICIENTS 
USED  BY  HUTCHINSON 


(D  ,D  ,D  )  = 


Et3 


-  (DXX’DVV’DVV>  *  <QvvQw>Qvv>  =  t(Q  ,Q  ,Q  > 


XX’  xy’  yy'  12(l_v2)  '  "xx  *  “xy  ’“yy '  ’  "<xx  ,v<xy  ,vyy '  ”  1  "<xx  >  vxy  ,yyy ' 

<Hxx’Hxy  ,Hyy5  “  Et  (fIxx  ,fIxy  ,fIyy 5  ’  ( Bxx ,Bxy ,Byx ,Byy 5  =  Et  (Bxx  ,Bxy  ,Byx  ,Ryy^ 


Bxx  =  asV1+“r)/oo 

Byy  =  arV1+aS)/ao 

13  s  va  a  y  /a 

xy  s  rTr  o 

Byx  s  '"Wb/bo 

Bxx  =  1+Bs  + 

+  [  12  ( 1-v2  )a  (1+a  )y2]/a 
s  r  s  o 

Byy  =  1+Br  + 

+  [12(l-v2)ar(l+as)Y2]/ao 

D  =  1  + 

xy 

+  [12(l-v2)vasarYsYr]/c*o 

^yx  =  ^xy 

Rxx  =  [lt"s(1-v2>’“0 

H  s  [ 1+a  ( l-v2 ) ] a 

yy  L  r  J  o 

Hxy  =  (l+v)-v/oo 

Byx  =  «xy 

Q  =  va  y  /a 

Xxx  s  s  o 

0  =  va  y  /a 

yy  r  r  o 

QXy  8  -CasYsCl+(l-^2)ar)  ♦ 

+  a_Y  ( l+(l-v2 )a  ) ]/2a 
r  r  so 

Qyx  +  Qxy 

as 

=  As/dst 

,  B_  =  EI  /Dd 

8  S  S 

’  Ys  =  V* 

ar 

-  Ar/drt 

.  Br  =  EIr/Ddr  , 

>  Yr  =  er/t 

a 

=  (1+a  )(l+a  )  -  v2a  a 
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Substituting  Eq.  (10)  and  (11)  into  (7)  and  (8),  leads  to: 

.3 


Et' 


12(1-0 


1  4p^?at2  2pQx 

27  ",xxxx  +  R  P!xx  +  - 003  R“  +  atW%x  =  0 


W* 


(13) 


— _J - -  -p*  _  L  w*  -  o 

Et ( l+ar )  ,xxxx  R  ,  xx 


Substituting  Eq.  (12)  into  (13)  yields  the  coefficients  A  and  B: 
12(l-v2)5ot5R2E(1+a  ) 

A  = - j - —  (14) 

O  p  p  p 

-  48(1-v  )  £ot  p  R 

B  =  - g - —  (15) 

where 


C  =  4E[4pit3  +  3(l+a  )R2(l-v2)t  - 


3ot(l-v2)  PqR2. 

E  ■ 


The  deformation  just  prior  to  buckling  is  governed  by 
Eq.  (11).  When  buckling  starts,  the  modal  deformation  pattern 
becomes  nonaxisymmetric .  The  modal  pattern  may  change  as  the 
postbuckling  process  continues.  The  initial  postbuckling  state 
is  assumed  to  be  of  the  form: 


W 

P 


vaR  2poX 

— g—  +  B  cos  — +  w 


,2  2p  x 

-  +  A  cos  + 


Et' 


where 

c 


3  ( 1-  v* 


) 


(16) 

(17) 
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and  w  and  f  are  infinitesimal  buckling  modal  quantities  repre¬ 
senting  the  displacement  and  stress  functions.  Substituting  Eqs. 
(16)  and  (17)  into  (7)  and  (8)  yields  the  following  compatibility 
(18)  and  equilibrium  (19)  equations  for  the  infinitesimal 
buckling  modal  quantities  w  and  f: 

.2  2t2(1+a  (1+v)) 

^  -p  -I-  ^  _p 

c'('l+or)  ,xxxx  4(1  ToT )  ,xxyy 

t2(l+ar(l-v2))  ^  tarYr 

+  c(i+or)  f,yyyy  +  O+aJj  w,xxyy 


vt  a  y 
r  r 

1  +  a 


w 


1 


-yyyy 


4p 

R  w,xx  +  ~^2  ^  5t+r) 


2P0X 

*  COS  (— R~)  W,yy  =  0 


(18) 


Et 


2Et 


Et' 


12(1-V2)  ’xxxx  12(1-V2)  W»xxyy  +  1  2  ( 1  -  v2  ) 


[1+6, 


+ 


12(1-v2)“rYr  1 
1+a -  ]  W 


a  y  _,4 
_  (  r  r )  Et  ~ 

»yyyy  Vi+v  c  ,xxyy 


+  vfe  ^  f 


1 +ar  c  ,yyyy  R  c  ,xx 


Et 


3  16p( 

1  2(l-v2) 


r  cos 


2p  x  Et(l+a  ) 

Jt-q  r 


TT 


T  COS 


2pox 


H 


4  PoCat2  2pox  4poE  5t4  2pox  „ 

+  ,2  cos  R  n2  c  C0S  R  f,yy 


(19) 


R 


R 


4p2rot  2p  x 

+  —~2 -  COS  -R-  +  atWjXX 
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cos 


+  ( 1  +  a  )  Etr  COS  — 5—  w  =  0 

r'  R  ,yy 


2p  x 
o 


where 


r 


48(1-v2) £ot2R2p2 


C 


Koiter  [ 1 1 ]  argues  that  for  the  unstiffened  shell  the 
nonsymraetric  buckling  pattern  can  have  the  form: 


N 

w  =  t  l  C  cos  [(2m-l)p  x/R]  cos  ny/R 
i  m  o 


(20) 


m=1 


He  suggests  that  the  prebuckling  stress  function  F  (Eq.  11)  is 
periodic  and  therefore  the  prebuckling  stresses  are  periodic, 
which  implies  a  periodic  distribution  of  circumferential  tensile, 
or  circumferential  compressive  membrane  stresses.  "Buckling  in  a 
nonsymmetric  mode  will  be  stimulated  in  regions  of  compressive 
circumferential  stresses  and  hampered  in  regions  of  tensile 
stresses.  We  may  therefore  expect  the  existence  of  an  asymmetric 
buckling  mode  with  nodal  lines  where  the  circumferential  stresses 
attain  their  maximum."  The  axial  period  of  such  a  mode  is  twice 
the  period  of  the  axisymmetric  prebuckling  state  response.  Thus, 
at 

x  =  ±  ^  =  (2  k+1  )  irR/2  pQ  ,  where  <  is  an  integer 

simply  supported  edges  are  assumed.  This  condition  is  incompat¬ 
ible  with  Eq.  (10),  but  Koiter  suggests  that  it  is  of  no  great 
consequence,  provided  that 
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Substituting  Eq.  (20)  into  (18)  results  in  a  fourth  order 
partial  differential  equation  in  f  with  constant  coefficients. 
The  solution  is  of  the  form: 


»  A. C  +  A„(C  .  +  C  , , ) 
^  j-  1  m  2r  m-1 _ m+1 

m=0  A3 


(2m-1)px 
]  cos  — 


R 


cos  ^ 


(21  ) 


where 


A.  = 


An  = 


A„  = 


arYrt2(2m-1 )p^n2  t2varYrn4  t(2m-1)2p2 

TT+^J"  r4  r3 


( 1 +  a„ )  R4 


0  2  2 
2p  n 

—J—  [5t+B] 

t2  (2m-1 )4p4  _  2(l+ar(l+v))t2  (2m-l)2p2n2 

c  ( 1  +  ar )  r4  +  C(1+ar) 


R 


1+°p(1’v  >  t2  n4 


l+o 


R 


Since  Eq.  (20)  cannot  exactly  satisfy  either  the  compati¬ 
bility  equations  or  the  equilibrium  equations,  we  assume  f  can  be 
made  to  satisfy  the  compatibility  equation  provided  that  the 
coefficients  of  Eq.  (20)  are  determined  in  some  mathematical 
least  square  sense.  This  is  accomplished  by  substituting  Eq. 

(21)  and  (20)  into  (19),  resulting  in  an  error  function  £: 


e(x,y)  =  e( A . , C  ,  cos 


( 2m- 1 ) PQX 


nr 


R 


cos 


2P03 

R 


cos  ^^)  (22) 


First,  the  trigometric  function  in  the  axial  direction  can  be 
condensed  into  a  single  function.  Second,  a  Galerkin  procedure 
will  be  applied  as  follows: 
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2  wR  j  Jt/2 

J-l/2 


:(x,y) 


(2q-1 )p  x 


cos 


R 


cos  dxdy  =0  (23) 


M 

R 

Q.  =  •  • 


Thus  the  infinite  set  is  formed: 


[ 


Et 


l  c. 


(2q-l)V 


12(l-v2)  **  ^  R4 

4  _4 


o  +  2Et4  n£  j.  c  (2q-1  )^Pq 


12(l-v2)  R2 


q  R2 


+ 


Et 


12(1—  )  R 


_  .  1  2  (  1  -  v  )  a  y 

2 \  r?  <1+6r  ♦  - — — >  l  0, 


Et4  n2  v  A1Cq  +  A2(Cq-1+0q+l )  (2q-1 )2p^ 
TTTT^  ~o  L  3 - -a-  a  * - ?; — - 


( 1  +  otr  )c  ^2 


R 


,  *Vl  n4  ,  A1Ca  + 

1 1  +ar )  c  ^  a 

Et5  y  A1Cq  +  A2^Cq-1+Cq  +  1  )  (2q-l)2p2 

L  A  ’  -  ^ 


Rc 


3 

,2  2 


R' 


-  at 


2  I  Cn  (2q'l)  P°  ]  [S4±2i  < 


(q-j)2p^i 


*  R2 


2  v  q-j ) 2P 


sin 


+ 


R 


(q+j-i ) 2p 


sin 


( q+d~i ) 2p^i  4p^trn2 


R 


-O-}]  -  [_o 


R 


4—  X  °q 


4p^(5t+r)n2  r  A,Cq  +  A2(0Q.1+Ca+1) 

A3 


R4c 


R  (  ttH-2  )  f  R 

•  ~r~^  {Tq^pt 


(q-j-1 J2p, 


sm 


(q-j-i )2Pql 

R - 
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r  _  (q+j-2)2p  L 

+  ( <1+3-2  J2po  sin  R 


R  ,  (q-j+l)2p  L 

+  lq-3+1 ) 2pQ  Sin  n 

R  (q+3)2p  L 

+  -  -  —  sin  - ^ }  =  0 


'(  <1+3  )2p. 


R 


(24) 


If  one  retains  r  equations  and  sets  C  =0  when  q  >  r,  a 
system  of  r  homogeneous  linear  equations  are  obtained.  The 

determinant  of  this  final  system  is  set  equal  to  zero.  With  the 
notation 


X  = 


K  = 


aRc 

Et 

2Rc 


and 


t  = 


n2t 

To 


the  characteristic  equation  can  be  written  in  the  form: 


H1  X5  +  H2X2  +  H^X  +  H4  =  0 


(25) 


where 


H1  =  1 


Hq  =  -  [ 


Q‘ 


2h: 


H 


+ 


2  ‘r  _2t,2 


6  C5t2(l+ar) 


8K2  Q2K2(l+a  )  4K2 


2  '  r 

A 


4K‘ 


] 


H,  +  H,  [-\  +  - 3. _  +  JS.  +  CST2(1+C1r>  , 

5  6  L,„ir4  «  „4/,  .  N  -...4  +  _ 4  +  j 


32K  Q2K*(l+ar)  64K 


32K 


Q^K2 
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H4  - 

n2  w4  c  2„2  22  4  / . 

2  Q,  H5  c?t  H5  c  5  t  (1+ar)H?  i 

6  51 2K6  52Q2K6(l+ar)  8Q2K4  8K2 

H5=  [-K2x2arTrC  +  J5-  T4varYrC  -  K2(l+ar)J2 


- 

4K4  +  4K2t2  +  T4(l+er  +  4ar^rc  ) 

1  ar 

4K4  +  4K2T2(l+ar(l+v))  +  T4(l+ar(l-v2)) 

H6  - 

16K4  +  (l+or) 

H7  - 

_L  +  _L 

p  -r  p 

Q2 

II 

CM  KA 

cy 

52K4  +  56K2x2(l+ar(l+v))  4-  T4(l+ar)l-v2) 

From 

Eq.  (25),  three  roots  (the  minimum  being  the  buckling 

load)  can  be  obtained.  The  computer  program  PVRCB  solves  Eq. 
(25)  for  the  roots  X. 

2.2.2  Nonsymmetric  Imperfection 

Th  above  procedure  can  be  repeated  for  the  case  of  an 
imperfection  geometry  of  the  form: 


W  = 
0 

-  €.j  t  cos  2K^x  +  J^t  cos  K^x  cos  Kgy  (26) 

where 


K,  ■ 

TiRmqo/2 

in 

IV) 

11 

TrRnqQ/2 

q0“(R/t)1/2  [12(1-v2)]1/4 
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and 

ra  and  n  are  integer  wave  numbers 


The  first  term  represents  an  axisymmetric  mode  of  amplitude 
?2t.  Hutchinson  [9,12]  assumed  that  the  buckling  mode  would  take 
the  form: 


W  =  cos  2K^x  +  cos  K^x  cos  K^y 


(27) 


and  derived  the  characteristic  equation: 


X  —  D2A  +  D .j  =  0 


(28) 


where 


D,  = 


(B1G1  -  B4&5)W2  -  (B3G5  +  B4G4)W^  -  B5G4¥^ 


D2  =  (B2G1  +  G2B1  +  B3G3)W2  +  B4G3  +  Gib5 
+  GpB,- 


B,  = 


K2  +  K2)2 


K 


4(K^  +  V?2)2 


2  2 

2cK1K^1 
(K2  +  K2)2 


K1 

B„  = 


B  3  =  -  c  K2  [  g 


2K, 


Oq  +  K2) 


+  il 


2  n 2  '  4- 


B  A  —  C  Kr)  5o  [ 


K, 


+  1  ] 


2  ^2L(K2+K2)2T4 


Br-  = 


K?«2 
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G,  = 


8Kt  +  i 


s0  m  k; 


G,  = 


-  4  KT?i 


G.  = 


-  c  K|[- 


K 


(Kf  +  K,) 


+  -1 
2x2  +  8J 


G[r  =  -  C  k£U 


K 


+  1] 


2  2(Kf  4-  K|)Z  '  4 


Hutchinson  has  shown  that  for  K1  =  K2  *  1 /2  a  minimum  A  is 
obtained.  Equation  (28)  is  also  solved  for  A  by  program  PVRCB. 
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2.3  Program  PVRCH 

This  program  calculates  the  Koiter  imperfection  parameter 
and  load  factor  for  a  ring  and  stringer  stiffened  cylinder  with 
isotropic  material.  The  stiffener  properties  are  assumed  to  be 
smeared  as  in  the  previous  case.  Out-of-plane  bending  stiffness 
of  each  stiffener  is  ignored,  but  the  torsional  rigidity  is 
included.  The  applied  load  can  be  either  an  axial  stress  or 
external  pressure.  Simply  supported  or  clamped  boundary  condi¬ 
tions  can  be  accommodated.  Only  an  axisymmetric  imperfection  is 
considered . 

Hutchinson  and  Frauenthal  [4]  analyze  the  imperfection 
sensitivity  of  an  axially  loaded  stiffened  cylinder.  The  general 
outline  of  the  procedure  has  been  covered  by  Budiansky  [13].  To 
arrive  at  a  solution,  Budiansky  states  that  three  basic  condi¬ 
tions  have  to  be  satisfied:  (1)  axisymmetric  prebuckling 
deformation  with  zero  imperfection,  (2)  nonaxisymmetric  bifur¬ 
cation  from  the  prebuckling  state  using  a  classic  linear 
bifurcation,  and  (3)  initial  postbuckling  behavior. 

Axisymmetric  Prebuckling  Deformation 

Equations  (7)  and  (8)  are  reduced  by  assuming  zero 
imperfection: 


lh[e]  -  LQ[W] 
ld[w]  +  LQ[F] 


0 


(29) 


FW  -FW  +2FW 
>xx  , yy  ,yy  ,xx  ,xy  ,xy 


0  (30) 


where  the  operators  are  defined  by  Eq.  (9)*  For  axially 
stiffened  shells,  the  Airy's  stress  function  is  related  to  the 
membrane  stress  resultant  within  the  shell  (Nx)  and  the  stringer 


N  +N 
x 


s 


F 

»yy 


p 

,xx 


p 

xy 


N 


X 


=  A  P  +  A  P 
xx  , xx  xy  ,yy 


+  B  W 
XX  ,  xx 


(31  ) 


The  evaluation  of  the  imperfection  sensitivity  of  a  cylinder 
subjected  to  an  axial  load  requires  the  determination  of  two 
parameters:  Pc,  the  critical  load,  and  b,  Koiter's  imperfection 

parameter.  The  expansion  of  the  critical  load  is  given  by  Eq. 

Q 

(1):  p/pC£  =  1  +  ea  +  e  b  where  PcJl  is  the  classical  load.  The 

imperfection  sensitivity  parameter  b  as  derived  in  Ref.  [4]  is 

b  =  f^2)  *(y(1),w(1))  +  2 

P  [  F^  *(W^  ^,W(1  J)  +  2  PU^  *(W°,WU))]  ^32^ 

where  the  notation 


A*(B, C)  =  /[A  B  C  +  A  B  C 

J  »xx  ,y  ,y  ,yy  ,x  ,x 

-  A  (B  C  +  B  C  )  ]  ds 


(33) 


1  xy 


y  ,x' 


It  is  assumed  that  in  the  neighborhood  of  the  bifurcation 
point,  the  deformation  and  stress  functions  can  be  expanded: 


W  =  W°  +  6W^  ^  +  62W^  +  .  . 


P  =  P°  +  SP^1 ^  +  62P^2^  +  . . . 


(34) 


where 


6  =  amplitude  of  the  imperfection 
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The  prebuckled  solution  is  expanded  about  the  critical  load,  P  , 
so  that 


w°  =  <  +  +  Z  (p-p0>2  »o  +  ••• 

p°  =  +  (P-P0)F°  +  f  (P-Pc)2  p°  +  ... 


(35) 


where 


C)  -  iLJ. 

V  '  3P 


P=P 


For  the  axisymraetric  prebuckling  state,  the  deformation 
pattern  is  assumed  to  be 

,  rO  O 

W  =  w 

and  the  stress  function 


F 


o 


p 

2ttR 


+  r 


(36) 


where  P  is  the  applied  load,  and  w°  and  f°  are  the  infinitesimal 
axisymmetric  modal  quantities  representing  the  displacement  and 
the  stress  function. 

For  the  axisymmetric  prebuckling  state,  the  governing 
equation  is  obtained  by  substituting  Eq.  (36)  into  (30),  which 
yields : 


^Dxx  +  H 


Q  iv  2Q 

xx^.o  ,  r  xx 

Jw  +  >-RH 


XX 


XX 


2ttR 


]wC 


w 


2 

R  n 


vA  P 
_ 


XX 


2ifEtH  R 

XX 


(37) 


and 


Hxxf 


O’ 


=  Q  w 
^xx 


0"  .  _  vAxyP 


+ 


R 


2nEtR 


(38) 


The  boundary  conditions  that  can  be  imposed  are  either 
simply  supported  or  clamped.  The  computer  program  uses  half  the 
cylinder  length  so  that  at  boundary  1  either  clamped  or  simply 
supported  boundary  conditions  are  imposed,  and  at  the  other  end 
symmetrical  conditions  are  assumed. 
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For  a  one-dimensional  system  of  equations  with  equal  nodal 
point  spacing,  the  derivatives  are  replaced  by  the  difference 
operators: 


and 


d  (  )  =  (  ^+1  (  *1-1 

dx  2  a 


where 


(  ) 


i+1 


2(  >1  *  (  >1-1 


2A 


A  =  L/N 

N  =  number  of  finite  difference  stations 


(39) 


The  resulting  system  of  algebraic  equations  is  solved  using  the 
standard  Potters'  method  [14]  of  forward  elimination  and  backward 
substitution.  Equation  (34)  can  be  solved  for  f°”  after  w°  has 
been  determined. 

Substituting  Eq.  (35)  into  (30)  results  in  an  equation 
similar  to  Eq.  (37): 


and 


rn  ^xx-,  .oiv 
[Dxx  +  — ]wo  + 

XX 


2Q 


XX 


■RH 


XX 


P  "  w° 

rc  .o  wc 

■3— S’Dw  +  -~— 

2  R  r2h 


vA 


xy 


XX 


2  irEtH  R‘ 

XX 


(40) 


2  irR 


Hxxfc 


=  Qxxw” 


w 


vA 


c  XV 

R  -  -  - 


2  irEtR 


(41) 


The  same  imposed  boundary  conditions  and  the  same  solution 
procedure  used  to  solve  Eq.  (34)  is  repeated. 
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Classical  Linear  Buckling  Problem 

Expansion  of  the  buckling  mode  (Eq.  34)  can  now  be  consi¬ 
dered.  For  the  linear  term  of  the  expansion,  the  governing 
equations  are: 


LD[¥(1 }]  +  LQ[F(1 }] 

lh[f(1)]  -  lq[w(1)] 

The  solution  of  Eq.  (42) 


2irR 


¥ 


(D 

,  xx 


-  r 


¥ 


(D 

,yy 


,0) 


w:  fv ' '  =  o 

c  ,yy 


+  w 


r(1) 


(42) 


¥v  7  =0 

,yy 


is  assumed  to  be  of  the  form: 


¥  ^  =  w  ^  cos  ny/R 

(D  .(D 


(43) 


F 


=  fK  cos  ny/R 


)  and  f ^ )  are  determined  using  the  previously  described 

o"  o" 

solution  procedure.  Note  that  f  and  w  must  also  be  obtained. 

c  c 

Initial  Postbuckling 

The  governing  equations  for  ¥^^,F^^  in  Eq.  (34)  are: 


LD[w(2)]  +  l8[p<2>]  =  jSg  W III  -  f°"  wj2>  -  »°"  ?l22 

(2ny/R) (f ^  ^  w^1 ^ 


=  ^  (H)2  [(f(l)w(D") 


2nR  ,xx  c  "  ,yy  "c  ,yy 

ti 


+  cos 


+  -  Ff^V1)’)] 


(44) 


,(2) 


LH[F^']  -  Ln[¥^']  +  ¥ 


(2) 


o"  . . ( 2 )  1  ,ns2  r 1  /  (1 )  (1 )" 


QL  J  c  ,yy  2  VR '  L2 
+  cos  (2ny/R)(w(l)"w(1)  -  w(l)'w(l)')] 


(Sr  c±  u"  v1'  ) 


(45) 
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The  solution  for  w(2)  and  p(2)  may  ^e  written  as 


W 


P 


(2) 

(2) 


=  wa  +  ¥g  cos  (2ny/R) 
=  fa  +  fp  COS  (2ny/R) 


(46) 


Prom  the  boundary  conditions  at  each  end  of  the  half  length 

shell  and  Eqs.  (44)  and  (45),  the  parameters  w  ,w„,f  and  f„  can 

a'  3’  a  3 

be  determined;  again  using  the  previously  described  solution 
procedure . 

Measure  of  Imperfection  Sensitivity 

In  order  to  determine  the  imperfection  sensitivity  of  the 
cylinder  due  to  an  imperfection  having  an  amplitude,  6,  the 
asymptotic  formula  given  by  Koiter  [ 1 5 ]  results  in  the  form 


=  1-5 


V 7 


ip  (s/t)2 


where 


_  b[P°*(w^\w +  F^1  ^(w°,w^  ^)]2 

(p°)2[P°*(w^  ;,wU  j)  +  2  PU  ;*(w°,w(1  ^)]2 
Thus,  when  the  state  of  stress  is  pure  membrane,  U  =  b. 


(47) 


(48) 
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2 . 4  Program  PVRCK 

This  program  calculates  the  buckling  load  factor  and  the 
Koiter  imperfection  sensitivity  parameter.  A  circular  cylin¬ 
drical  layered  orthotropic  shell  is  assumed.  Both  ring  and 
stringer  stiffeners  are  smeared  as  in  the  previous  case.  Out-of- 
plane  bending  stiffness  is  ignored  but  the  torsional  rigidity  is 
included.  The  applied  load  can  be  any  combination  of  axial 
stress,  external  pressure,  and  torsion.  The  boundary  conditions 
are  not  specifically  imposed,  but  simply  supported  conditions  are 
approximated  by  the  selection  of  the  radial  displacement  response 
functions . 

Khot  has  studied  the  buckling  and  postbuckling  behavior  of 
composite  cylindrical  shells  [5,16].  He  also  has  investigated 
the  imperfection  sensitivity  of  these  shells.  Several  computer 
programs  were  developed  by  Khot.  These  programs  have  been 
revised  to  more  efficiently  evaluate  the  buckling  characteristics 
under  either  separate  or  combined  states  of  stress  resulting  from 
axial  compression,  external  pressure,  or  torsion.  Some  develop¬ 
ment  of  the  theory  has  already  been  published,  but  many  key 
points  remain  in  a  preliminary  and  unpublished  form.  Some  of 
these  data  have  been  graciously  supplied  by  Khot. 

The  shallow  shell  strain-displacement  equations  and 
Donnell's  equilibrium  equations  (Eqs.  (2)  and  (3))  are  again 
used . 

The  constitutive  equations  take  the  usual  form: 


(49) 
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e 

K 

X 

X 

e 

i  y 

*  ;  <  =  < 

i 

KJ  > 

I  e 

! 

1  K  ! 

xy 

!  xy 

1 

J 

l  J 

C*m  =  I  BL<hi+1-hi) 

KHid  =  7  l  Btm  (hl+1~hl) 

D*n,=5  l  B*«.(hi+rhl’ 

=  elastic  moduli  of  each  layer 
hi  =  distance  to  each  layer  from  reference  surface 

Note:  A^m,  K^,  and  can  he  augmented  by  the  smearing  of  the 

stiffener  properties. 


Solving  for  e  in  the  first  equation  of  (49)  and  substituting 
into  the  second  results  in  the  semi— inverted  form: 


d*  =  D  -  KTaK 


(50) 


Using  the  definition  of  the  stress  function  from  Eq.  (6),  we  can 
obtain  the  equilibrium  equation  in  the  form: 
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d,_  F  +  (2d,  -d. , )F  +  (d . . +d__-2d,, )F 

12  ,xxxx  62  16  , xxxy  11  22  66  ,  xxyy 

+  (2d,,-d0,)F  +  d„,F  +  F  (w  +  A)  (51) 

61  26  , xyyy  21  ,yyyy  ,xxv  ,yy  R'  w 

-  2F  w  -  d*  w  -  4d.,w  -  (2d*  +4d*,)w 

, xy  ,xy  11  ,xxxx  16  ,xxxy  v  12  66'  ,xxyy 


-  4d*,w  -  d*  w  -  p  =  0 

26  , xyyy  22  ,yyyy  * 


The  compatibility  condition  is  written  as: 


a22P,xxxx  2a26F,xxxy  +  ^2a1 2+a66^F,xxyy  2al6F,xyyy 
+  a1 1P*yyyy  ~  ~  d12w,xxxx  "  ^2d6l _d26^w,xyyy 


w 


(52) 


d21w,yyyy  +  w,xy  w,xxw,yy 


,  XX 


R 


-  (2d,  -d  ,)w  -  (d  +d0_-2d,,)w 

62  16  ,  xxxy  11  22  66  ,  xxyy 


For  the  above  two  equations,  the  subscript  notation  of  1,2  and  6 
correspond  to  the  material  properties  in  the  two  orthogonal 
directions  and  shear,  respectively. 

After  a  lengthy  normalization  procedure,  Eq.  (51)  can  be 
written  as: 


kF  +  vF  +  iFF  +  £F  + 

, xxxx  , xxxy  Y  , xxyy  ,xyyy 


+  F  w 

+  i  F 

—  2F  w  + 

,xx  ,yy 

2  ,xx 

» xy  , xy 

+  A  w 

+  £W 

+  w  - 

2  , xxyy 

»xyyy 

4  ,yyyy 

and  Eq.  (52)  can  be  written  as: 


XF 

,yyyy 

4-  w  + 

4  ,xxxx 
P<t>  =  0 


vw 

,  xxxy 


(53) 


4 

<!> 


F 


+  $F  +  4aF  —  yF 

xxxx  ,xxxy  ,  xxyy  '  , 


+ 

KW 

+  VW 

,  XXXX 

,  xxxy 

+ 

2w,xx  = 

w^xy  -  « 

+  tyw  + 

,  xxyy 


,  xxw ,yy 


+  4d)F 

xyyy  ,yyyy 

£w  +  Tw 

>xyyy  ,yyyy 


(54) 
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where 


d1 2  (al1,l/2  1 

K  *  a22  ^*22  * 


5  = 


- m  <^)1/4<*),/2 

(a00d*  )1/2  a11 


2d6l"d26 


22  22 
a 


Y  =  8 


( a  1 1  a2  2 ) 


12  m1/2 

i74  U) 


6  = 


o  26  m-i/2 

77I777  u> 


a  = 


2a1  2+a66 


c  = 


v  = 


*  = 


♦  = 


X  = 


(a,, 

a22 

jm 

d?6 

d* 

q22 

)1/4 

d* 

26 

’  (dn 

d*^ 

a22 

jm 

(2d 

62" 

d,6> 

( a22 

d22 

)1/2 

d11+d22 

-2d66 

( a22 

d22 

yr/2- 

(a11 

a22 

dTi 

d22 

1/2 

) 

d21 

(a1 1 

d22 

772 

»i  1 . 1 /4 


(^) 


‘  a 


(♦) 


-1/2 


22 


Eqs.  (53)  and  (54)  are  now  of  the  form  that  the  critical  load 
parameter,  A  ,  and  the  imperfection  sensitivity  parameter  of  the 
shell  can  be  derived. 


31 


For  the  present  case  a  solution  can  he  obtained  with 
Eq.  (34)  rewritten  as: 


F  =  AcF°  +  6F^ 1 ^  +  52F^2^  +  **• 


w  = 


+  5V2>  +  .. 


Substitution  of  these  equations  into  (53)  and  (54)  and 
collection  of  like  powers  of  6  results  in: 


,(D 


(D 


Li  [wV  '  ']  -  LjF(l  }]  =  A  [F 


w 


(1) 


+  F°  W(1 l  -  2P°  wll!]  (55a) 


c  , xx  ,  yy 


(1) 


>yy  ,xx 


*xy  ,xy- 


,(2) 


l2[wv  ']  +  l5[f^;]  =  0.0 
and 


(55b) 


Lll>(2)]  -  L9[F^]  -  A  [ F° 


(2)  _,o  (2) 

MLW  J  -  i/oLr  J  —  a  |  r  W  +  P  W 

1  2 L  cL  ,xx  ,yy  ,yy  ,xx 


2F°  w(2)] 

.xy  xy  J 


=  F^w^  +  F^w^  -  2F*'1'^w^1^ 

»XX  , yy  ,yy  ,xx  ,xyw,xy 

+  1,  [f^2>]  =  (w(lb2  -  J'U') 


iXy ' 


w  ,  '  w 
»xx  ,yy 


(55c) 


(55d) 


where 


y  ] 

L2[  ] 


y  ] 


=  ~r  i 

4l  j 

=  <[  ] 


$ ] 


,  xxxx 


+  n[  ] 


,  xxxx 

+  X[  ] 


+  #[  ] 

,  xxxy  2L  J 


+  v[  ] 


xxxy 


+  *[  ] 


+  2[  ] 

>yyyy  ,xx 


+  rT  I  +i[  1 

xxyy  ■’L  J,xyyy  4L  J,yyyy 

+  5[  ] 

, xxyy  , xyyy 


—  3  [  ]  +4ot[ 

, xxxx  L  J , xxxy  L  J , xxyy 

-  y[  ]  +  4  4> L  ] 

>xyyy  ,yyyy 


Neglecting  the  nonlinear  coupling,  Ac  can  be  determined  for  the 
case  of  combined  axial  compression,  external  pressure  and 
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torsion.  It  is  assumed  that  the  classical  radial  displacement 
under  torsion  can  be  of  the  form: 


w  =  h  sin  SI*  cos  £  (y-xx)  (56) 

where 

(a22d?2>1/2 

r  =  R 

(4R2d*2a22)1/4 

(4R2d|ia22)1/4 

x  is  an  integer 


If  the  notation 


M  =  ^  “  x 

l  r 


N  =  - 
r 


T,  mir  n 
p  -  —  -  ?  T 


is  used,  then  the  classical  radial  deflection  for  the  initial 
buckling  mode  can  be  written: 


(1  )  h  . 

T  =  —  Cl 


w'-'  =  sin  (Mx-Ny)  +  sin  (Px  +  Ny) 


(57) 


or 


w  ^ 1  ^  =|  (a1+a2) 


It  is  assumed  that  the  Airy’s  stress  function  can  be  expressed  in 
a  similar  form: 


)  =  P  A  +  F  A 
11  2  2 


(58) 
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Prom  the  compatibility  equation  (54)  and  the  definition 
of  w(1)  and  we  obtain: 


Pi  - 


h  kM^  -  vM3N  +  ^M2N2  -  SMN3  +  XN^  -  2M2 

If  4  IJ  3 — - T~ 2 - 3 - n“ 

-M  +  3M  N  +  i)oM  N  +  yMN  +  4<f>N 


T 

h  _3 

2  Tc 
0 


(59a) 


4  3  22  2  —  li  o 

p  h  kP  +  vP  N  +  U>P  N  +  £PN°  +  Xn  -  2P^ 

2  2  4  4  3  ,  2  2  3 - IJ — 

-  P  -  3P  N  +  4aP  N  -  yPNJ  +  4 <f>N 

T 

=  -  h  _4 

2  T«r 


(59b) 


Prom  Eq.  (53)  we  obtain: 


T1A1  +  T2A2  "  T3P1A1  "  T4P2A2  =  "  Xc  I  ^Nx(m2ai  +  p2A2  > 


+  N°N2(A1+A2)  +  2  N°  (-NMA1  +  NPA2 )  ] 


(60) 


where 


T 


1  =  |  (-$  -  nP3N  +  |  M2N2  -  ?MN3  +  4  N4) 


T2  =  I  (4  pZ<  +  np3w  +  §  p2fj2  +  4PN3  +  |  P4) 


Collecting  like  terms,  we  can  now  solve  the  critical  load 
factor: 


2  2 
To  T„ 

-  r  m  +  T  +  — —  +  — — 1 
L  1  l2  T*  TVJ 
>  =  _ 5  o 

c  „o  /m.2.„2 


N°  (M^+P^)  +  2  N°N2  +  2  N°  N(P-M) 
A  y  xy 


(61) 
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where : 


(a1 1a22>1  ^ 

4(d|2)l/2 


R 


are  the  stress  resultants  calculated  from  the  applied  loads 
Nx,Ny  and  Nxy  and  are  assumed  ^ 0  a  pure  membrane  state. 


In  order  to  evaluate  the  coefficients  a  and  b  of  the  expan¬ 
sion  of  Eq.  (l),  Eq.  32  must  be  used  with  Eqs.  (57),  (59a),  and 
(59b).  The  coefficient,  a,  has  been  determined  to  be  zero  from 
consideration  of  the  periodicity  of  the  circumferential 
displacement . 

With  F^  )  and  w^  )  determined,  the  postbuckling  w(^)  and 

(?) 

Fv  '  are  assumed  of  the  form: 


•*  ■*■  ^  *  * 


1  .  2ny 

+  7T  sm  — *- 
d  r 


1 


P(2)  =  l  B  si„  +  i 

i=odd  1  * 


i=odd 
1_ 
2 


1  .  2ny 

+  2  sin  r 


where 


l 

i=odd 


M.  =  il  + 

1  z  r 


P  =  ll  _ 

i  i  r 


cos  £ 

i=odd 

Y • ( cos  P . x  - 
1  1 

cos  l 

1  i=odd 

6 . ( COS  P . X  - 
1  1 


Yi(sin  M^x  +  sin 

cos  M.x) 
x  ' 

<5^(sin  M.x  +  sin 
cos  ]VLx) 


PjX) 

(62) 
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It  remains  to  determine  the  coefficients  a.,  3.,  y.  and  6  .  To 

111  i 

obtain  relationships  for  these  coefficients,  Eqs .  (59a),  (59b) 
and  (62)  are  substituted  into  Eqs .  (58)  and  (55c)  to  obtain: 

4L  l  “i  sin(1T£)  +  ^  CSR  l  T;iA4  +  -Tj-  SNR  l  y^J 

+  h[l(fi)  SNR  l  y.^  +  1  CSR  l  y.A3] 

+  CSR  l  TiA2  -  i(|^)2  SNR  l  y±n2] 

+  c[f(jp)3  SNR  l  Y.fi1  +  1(|^)3  CSR  l  y.Aj 

+  CSR  l  y.Ao  +  l(g^)4  SNR  l  y.flj  (63) 

-  KtlPiCj1)4  sin(-^p)  +  1  CSR)  l  5iA4  +  1  SNR  l  6^] 

-  ^(f1)  SNR  l  6.^  +  1(|H)  CSR  l  6.  a3] 

+  ^(f1)2  CSR  I  6.A2  -  1(§^)2  SNR  I  6.n2] 

-  SNR  I  «i81  -  ^(f1)3  CSR  I  6iA1] 

-  X[l(|^)4  CSR  I  6.Ao  4-  1(§^)4  SNR  I  6.flo] 

-  2[  J  a.(-^)2  +  ^  CSR  I  6.A2  +  1  SNR  l  6.^] 

"  Xc[Ny{4^2  CSR  l  Vp  -  ^F1)2  SNR  I  W 


+  N^  {-  l  (-£— ) 2  sin 


17TX 


-  CSR  l  y±A2  +  1  SNR  I  y^} 

+  2Nxy{-  Z<|r>  SNR  l  Vi 

+  |  (|£)  CSR  I  yiA1  }] 

=  -  j-l?1  +  n2[m+p]2  A  k2 

5  6 
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where 


SNR  =  sin(~£) 

CSR  =  cos(^) 

Aq  =  sin  (M^x)  +  sin  (P^x) 

sin  (M^x)  +  sin  (Pix) 

Ap  =  M?  sin  (M.x)  +  P?  sin  (P.x) 

A^  =  -  M?  sin  (M^)  +  P?  sin  (Pix) 

A^  =  Ml  sin  (M^)  +  pf  sin  (P^^x) 

«  =  cos  (M.x)  +  cos  (M.x) 

o  1  1 

=  -  M.  cos  (M.x)  +  P.  cos  (P.x) 

1  l  l  i  i  7 

^2  =  ^i  003  ~  P2  cos  (Pjjx) 

=  M^  cos  (Mix)  +  P?  cos  (P.x) 

^  =  -  M4  cos  (M^x)  +  p|  cos  (Pix) 

Also,  Eqs.  (59a),  (59b)  and  (62)  are  substituted  into  Eq.  (55d) 
to  obtain: 

k[I  “.(ll)4  sin  (ip)  +  1  CSR  I  y.A4  +  1  l  y.Q4] 

+  v[l(|£)  SNR  I  Y^  +  |  (§2)  CSR  £  Y±  A^ ] 

+  ^(f1)2  CSR  l  YiA2  -  l(fl)2  SNR  l  Yin2]  (64) 

+  SNR  l  Yi01  -  i(§^)3  CSR  l  y.A^ 

+  )4  CSR  I  y.  Aq  +  l(fl)4  l  y.%] 

+  2[-  l  a.(li)2sin  (ip)  -  i  CSR  I  YiA2  +  i  SNR  I  Y-^] 
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+  f[I  6t  (i^)4  sin  (i|i)  +  1  CSR  I  S1A4  +  1  SNR  l  J.a^ 

-  6[i(|p)  SHR  I  «.a3  +  1(|H)  CSR  l  S.Aj] 

+4o[j(|2)2  CSR  l  SjAg  -  jig2)2  SHR  J  5i!)2] 

-  y[j  (|2)5  SHR  J  qa,  -  1  (|^)3  CSR  I  S1A1  ] 

+4*[|  (|S)  CSR  I  5lA0  +  i  (|S)4  SNR  J  S.aJ 
2 

=  -  £-  n2[m+p]2  a^2 


Multipling  Eq.  (63)  by  sin  and  integrating  between 
we  obtain: 


2  irR 


|L«i(x;)‘t]  -  <L 0^ <x")4  -  201(ir>‘i]  +  ‘c<“i(r> 

2i 


■  i  TT  \  4 


■  i  ir  s  4 


’  i  TT  \  2 


i  TT  \  2 


-  h2lh  ,  T4)fI2  4*m2 

4  V  VN  «2  TI^?r 


(65) 


Multipling  Eq.  (64)  by  sin  and  integrating,  we  obtain: 


-i^)4  -  2«i(f)2  ♦  f  M^>4 


=  h2N2  m  77 


2i 


.2  /  •  2  :  2s 

l  ( 1  -4m  ) 


Eliminating  g  we  obtain: 


a . 
l 


2h2N2m2  jri 
(i2-4m2)  £2 


{(^  +  *4)  4  +  *  -  2(^)2} 


T 

i 

T, 


ITT  ‘ 


(i1)4  +  ?  'i1)2  +  -  2>2 


and  from  Eq.  (66) 


,  2W2  2 
■h  N  m 


TT1 


=  4>[  o  0  - p  vq— 

(i2-4m2)  2  A2  1 


(:M4 


r  <  1  \2-,  1 

i {4  '  2  (T7} 


(66) 


(67) 


(68) 
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Multiplying  Eq.  (63 )  "by 


cos 


[sin  M.x  +  sin  P.x] 


1  .  2ny 

+  2  Sin  R 


[cos  P.x  -  cos  M.x] 

u  1  1  J 


and  integrating,  results  in: 


’ri[T5(Mi+Pi)  -  + 

+  Iff2)4]  -  +  |(|I1)2(Mrpi) 


nr2nWM3  ^3 


p/2h\ 2/„2 , ^2 


5/2nx3, 


|(^)3(Mi-Pi)  +  Itf1)4-  ^(M2+P2)]  +  xc  Ii[N°(|^)2 


7VR 


(69) 


+  ^(m2+P2)  _  NOy(2n)(Mi_Pi)]  .  £{b  + 


Multiplying  Eq.  (64)  by  the  same  factor  and  integrating,  results 
in : 


Yi[f(Mi+pi)  -  ?(fr)(M3-pi3)  +  |<l!1>2(4fPi)  -  f<ir>3(Mi-pi) 

„4,„4n  .  B/„3  „3- 


+  l^)4  .  ^  (M2+p2}]  +  5.[l(HpPp  +  |0«J-P’) 

+  2<t»(|^)4] 


(70) 


+  a(f^)2(M2+P2)  +  £(f£)(M?-pJ)  +  ^(fp)3(M.-P..) 


4  R 
2 


i  i 

2 


hi  „2  (M+P) 

4  i  TT 


Solving  for  6^  in  Eq.  (70)  and  substituting  into  Eq.  (69),  we 
obtain : 


Yi  = 


-4h2N2m2TT2  /I  \  f /T3  ,  ^4 

*2  (i^UT5  +  T6 


) 


T9  +  T8} 


T7T9  +  T8p 


(71) 
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,  r4h2N2m2Tt2  /Is  ,  1 

5-  =  -  ‘s -  <Ti>  +  V8]  17 


4  -  -  t; 


(72) 


where 


T  = 
7 


|(m4+p|)  _  n(|lL)(M3_p3)  +  §(§^)2(M2+P2) 

-  C  (~)3(M.-P.)  +  |(|^)4  +  ac[nJ{m2+p2} 


+  2N°  (§^)2  -  2N°y{M.-P.}  |S] 
k(M4+P4)  -  v(|S)(M^-P^)  +  (M2+P2) 

-  5(fi)3  (M.-P.)  +  2  A  (f^)4  -  2  (M2-P2) 

■9  ~  *  (Mi+Pl)  +  e(|I1)(M3-P3)  +  4«  (|^)2  (M2+P2) 
+  Y  (|^)3  (M.-P.)  +  8<p  (|^)4 


2n- 


T 

o  ~ 


T/,  = 


Equation  (32)  can  he  rewritten  as: 

b  =  j fi~F^ 

-  +w|;)W(2)}]axay 

+  F^)(W0))2  .  2P(2)w(1l)wO)]dx445'  / 

-  X 


'XX  ,y  ,y  *  ,yyw,x  w,x 


[F(2)(W(! ))2 


,  xx 


Kxx(Wly>)2  +  p°yy(»(i))2 


,  X 


,xyw!))k,!)W 


(73) 


Substituting  w( 1 ) ,  pO),  w(2),  and  p(2)  into  Eq.  (73)  results  in 
the  following  expressions: 
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a 


where 


T 


,w2  2  3* 

4N  m  u  10 

2 —  f — 

V  T11 


T  T 

1 0  =  +  t7  ^  ^  £  ai  ~2 — ? 

,u  x5  i6  i=odd  1  i-4m2 


-  Z  *i 


2i 


.  2  .  2 
l  -4m 


6. 

I  — 

L  l 


Y- 


I  f  1 

i=odd  1 


T 


11 


No{mi)2  +  n2t2}  +  non2  _  2NayN2x 


(74) 


With  an  imperfection  amplitude  in  a  given  harmonic  6  specified, 
the  knockdown  factor  ( A  / A  )  can  then  be  determined  by 

D  U 


(1  -  W3/2  -  i  <»./»„>  ^ 


(75) 
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2. 5  Program  PVRCA 

In  the  previous  sections,  the  computer  programs  are  rather 
limited  in  scope  with  respect  to  the  assumed  imperfection  modes. 
At  most,  just  two  imperfection  modes  could  he  considered  simul¬ 
taneously.  In  order  to  develop  a  correlation  to  experimental 
data,  general  imperfection  data  must  he  derived  from  specific 
measured  geometric  data.  Arbocz  has  written  a  procedure  that 
accounts  for  the  effects  of  a  general  imperfection  when  a 
cylinder  stiffened  hy  rings  or  hy  stringers  is  subjected  to  an 
axial  load.  Further,  the  effects  of  nonlinear  prebuckling  and 
boundary  conditions  are  incorporated  into  the  procedure. 

The  in-plane  and  out-of-plane  bending  of  the  stiffeners  are 
considered  in  the  analysis.  A  general  imperfection  geometry  can 
be  specified.  Either  load  or  end-shortening  increments  can  be 
considered . 

A  detailed  description  of  the  development  of  governing  equa¬ 
tions  for  the  multimode  analysis  is  given  in  Ref.  [l7]«  Arbocz 
and  Babcock  [6]  give  a  synopsis  of  the  procedure  and  show  how  the 
quasilinearization  method  of  Newton  is  used  to  solve  the  Donnell 
equations  (Eq.  7  and  8).  The  solution  to  these  equations  are 
based  upon  the  approximation: 


m+1 

=  wra 

m 

+ 

6  Wm 

m 

l 

m+1 

It 

B 

+ 

6  F 

m 

where 

¥m,Fm  =  mth  approximation  to  the  solution 


<$Wm,<5Fm  =  correction  to  the  mth  approximation 


(76) 


Substituting  into  Eqs .  (7)  and  (8)  and  neglecting  squared  terms 
yields : 
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(77) 


LH(6Fm>  -  V6V  +  6Wm,xx/R  +  %/VVV  = 

LQ(6V  +  V*V  -  5Fm,xx/R  -  LNl/V*V 


-  E 


(D 


m 


-  Lwt(W  +W  5F  )  =  -  E 
NL '  m  ,  m '  m 


(2) 


where 


Lm(S,T) 


(1) 

m 

(2) 

'm 


ST  -2ST  +ST 
.xx  ,yy  ,xy  ,xy  ,yy  ,xy 

W  -  W  +  Vxx/R  +  ^(Wm.V2W) 

LQ(Pm)  +  W  -  Pm’xx/R  -  WW"* 


The  radial  imperfection  is  assumed  to  be  represented  by: 


N, 


N, 


1  .  _  _ 
w  =  t  l  wirt  cos  ix  +  t  l  l  w 

k,£=1 

N„ 


1  io - —  ■  -  ,  4  l  "kA  sin  kx  003  *7 


(78) 


+  t  l  l  w 
k,  4=1 


k£  sin  kx  sin  4y  =  t  I<wLwJ 


where 


x  =-  ttx/L 
y  =  y/R  =  0 


Similar  solution  forms  can  be  made  for  the  four  unknowns  of 
Eqs .  (76)  and  (77). 


w 


m  =  - 


,  ff 

tv  XX 


1  +a 


4-  t  L  [  w  ] 
wL  J 


6w  =  t  L  [iw] 
m  wL  J 
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TABLE  3.  RING-STIFFENED  SHELLS 


Case  Number 

Parameter 

1A 

12 

,  21 

24 

R 

3.929 

3.929 

3.891 

3.891 

H 

0.0165 

0.0173 

0.0185 

0.0223 

l 

5.00 

11.00 

11.125 

8.90 

No  Rings 

9 

22 

22 

18 

d 

r 

0.50 

0.50 

0.50 

0 . 50 

e 

r 

-0.0193 

0.0221 

0.0226 

0.0241 

A  x  102 
r 

0.175 

0.391 

0.3072 

0.295 

I  x  106 
r 

0.0877 

0.271 

0.2156 

0.1961 

J  x  106 
r 

0.807 

0.90 

0.69 

0.62 

6 

0.0 

0.0 

0.00187 

0.00195 

P  Test 

cr 

486 

538 

480 

771 

ncr 

8 

8 

8 

8 

PVRCB 

470.3 

562.6 

531.2 

705.7 

ncr 

7 

10 

11 

11 

PVRCH 

400.5 

505 

544.2 

652 . 4 

ncr 

7 

9 

12 

11 

PVRCK 

469 

508 

531.7 

776 . 5 

ncr 

12 

9 

10 

10 
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TABLE  4 


STRINGER-STIFFENED  SHELLS 


Case  Number 

Parameter 

101 

105 

121B 

124A 

R 

3.694 

3.703 

3.759 

3.762 

H 

0.0349 

0.0194 

0.0227 

0.0180 

l 

10.89 

10.89 

5.07 

2.12 

No.  Stringers 

36 

36 

36 

36 

d 

s 

0 . 649 

0.649 

0.666 

0 . 666 

es 

0.0787 

0.071 

0.0521 

0.0498 

A  x  102 

s 

2.60 

2 . 60 

12.13 

12.13 

I  x  106 

S 

46.33 

46.33 

6.73 

6.73 

J  x  106 

s 

78.95 

78.95 

17 . 67 

17 . 67 

6 

0.0 

0.0 

0.0029 

0.00182 

P  Test 

cr 

7026 

3980 

2640 

3360 

ncr 

- 

- 

6 

7 

PVRCH 

5370 

2664 

3139 

7371 

n 

cr 

7 

7 

6 

6 

PVRCK 

4700 

2862 

2120 

2867 

ncr 

6 

4 

6 

6 
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the  critical  buckling  load.  If  too  small  a  step  is  taken,  no 
solution  may  result.  If  too  large  a  load  step  Is  taken,  a  root 
may  be  missed. 

Programs  PVRCB,  PVRCH  and  PVRCK  all  give  reasonably  good 
results  for  a  ring-stiffened  cylinder.  The  maximum  spread 
between  the  experimental  results  and  the  three  programs '  pre¬ 
dictions  is  18  percent.  For  the  stringer  stiffened  shells,  a 
much  greater  deviation  is  observed  between  the  reported 
experimental  results  and  the  two  programs  PVRCH  and  PVRCK.  It 
appears  that  the  effect  of  boundary  conditions  is  more  critical 
for  this  class  of  shell  than  the  ring-stiffened  shell,  as  would 
be  anticipated.  PVRCK  appears  to  be  the  more  conservative  and 
consistent  program. 

As  a  demonstration  of  the  output  from  each  program.  Case  24, 
a  ring-stiffened  cylinder,  is  used.  In  Vol.  II  (User's  Manual), 
the  CRT  displays  are  given  for  Programs  PVRCB,  PVRCH,  and  PVRCK, 
respectively.  The  results  from  PVRCB  require  no  further 
explanation.  From  PVRCH,  a  minimum  load  of  862.2  lbs  at  a 
response  mode  of  nCP  =  11  is  obtained.  In  order  to  obtain  the 
minimum  load  we  must  use  the  relation: 


1/3/  .00195^/3 


=  862. 2(1-1. 89(. 2791) 


=  652.4  lbs. 


From  PVRCK  a  minimum  stress  resultant  of  31.76  lb/in  was  pre¬ 
dicted.  The  minimum  axial  load  is 

P  =  N  2irR 


cr  x 


=  31.76(2)tt(3.891) 
=  776.5  lbs. 


The  experimental  critical  buckling  load  was  771  lbs. 
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The  program  PVRCA  cannot  be  run,  at  present,  on  an  inter¬ 
active  basis.  Therefore,  no  CRT  output  is  possible  with  the 
present  machine  configuration.  If  larger  core  limits  were 
available,  interactive  sessions  could  easily  be  obtained.  This 
program  uses  imperfection  data  from  existing  shells  and  then 
determines  the  critical  buckling  load.  Without  prior  knowledge 
of  what  response  modes  are  important,  the  selection  of  the  major 
responding  modes  is  by  pure  chance.  For  example,  the  data  for 
the  first  ring-stiffened  cylinder  ( 1 A )  was  used  with 
imperfections  of  the  first  two  terms  in  Eq.  (78)  being  0.00001 
for  modes  ( 9 ,0) , ( 1 8, 0) ,  (27,0) , (9,7) , (9,8) , (9,10) , (9, 1 1 )  and 
(9,12).  The  estimated  critical  buckling  load  is  X  =  1.234  or 
523  lbs.  For  Case  24,  where  the  imperfection  is  known  to 
exist  ( 5  =  0.00195),  the  predicted  critical  buckling  load  is 
774  lbs  with  the  same  responding  modes. 

A  further  demonstration  of  the  adequacy  of  the  program  PVRCA 
is  represented  by  an  example  of  the  third  category.  From 
Arbocz's  thesis  [ 1 8 ] ,  the  description  of  an  unstiffened  imperfect 
cylindrical  shell  is  given  (see  Tables  5  and  6).  The  form  of  the 
imperfection  is: 


N  N 

w(x,y)  =  l  l  A  cos  m  ^  cos 

m=0  n=0  mn  R  L 

N  N 

+  l  l  Bmn  sin  m  ^  cos 

m=1  n=0  mn  R  L 

C  cos  m  ^  sin  ^nrrx 
mn  R  L 

D  sin  m  ^  sin  2-?—x 

mn  R  L 

where 


N  N 

+  l  l 

m=0  n=1 
N  N 

+  I  l 

m=1  n=1 
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TABLE  5. 


A7  Fourier  Cosine  Coefficients 


\  M 

n\ 

0 

1/2 

1 

3/2 

2 

5/2 

3 

7/2 

0 

0.068 

0.068 

1 

2 

-0.402 

-0.457 

-0.518 

-0.049 

-0.155 

-0.085 

-0.103 

-0.054 

3 

-0.456 

-0.601 

-0.586 

-0.045 

-0.099 

-0.062 

4 

0.103 

0.048 

0.136 

5 

0.193 

0.073 

0.242 

0.056 

6 

0.114 

0.145 

7 

-0.047 

8 

-0.050 

9 

0.051 

10 

-  -0.048 

-0.073 

TABLE  6. 


A7  Fourier  Sine  Coefficients 


B. 


mn 


D 


mn. 


'''v  M 
Nx. 

0 

1/2 

1 

3/2 

2 

5/2 

3 

7/2 

0 

1 

2 

-0.241 

-0.348 

-0.308 

-0.271 

-0.100 

-0.126 

-0.054 

-0.079 

3 

0.150 

0.189 

0.066 

4 

0.504 

0.129 

0.644 

0.100 

0.200 

0.048 

0.112 

5 

0.128 

0.167 

0.047 

6 

-0.165 

-0.149 

-0.200 

-0.108 

-0.063 

-0.047 

7 

-0.074 

-0.075 

-0.056 

8 

0 . 069 

0.049 

9 

13 

-0.048 

& 


L 


mn 


B 


mn  ttL  i 


Q  _  _ £ 

mn  it  LI 


D 


mn 


4 


2irR 

/  0  Jo 

I  f2.R 

'o  Jo 

L  f  2uR 

o  Jo 

L  f  2irR 

'  o  Jo 


w(x,y)  cos  m  £  cos  — ^  dy  dx 


w(x,y)  sin  m  J  cos  ^nirx  dy  dx 


w(x,y)  cos  m  J  sin  dy  dx 


w(x,y)  sin  m  J  sin 


E —  dy  d* 


A  stiffened  cylindrical  shell  (A7)  was  used  with  the  follow¬ 
ing  properties: 

Radius  =  4-003  in. 

Length  =  8.00  in. 

Wall  Thickness  =  O.OO4494  in. 

E  =  15.1  x  10^  psi. 

v  =0.3 

The  largest  19  coefficients  were  used  to  predict  the 
buckling  load  and  was  found  to  be  X  =  0.6.  The  experiment  gave 
X  =  0.554. 
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4.0  BEHAVIOR  OP  KOITER'S  PARAMETER 

For  most  solution  procedures  used  in  the  previous  section, 
the  imperfection  sensitivity  of  a  shell  is  based  upon  the  evalua¬ 
tion  of  Koiter's  parameter,  b,  in  accordance  with  Eq.  (73).  In 
order  to  view  this  concept  in  proper  perspective,  it  was  decided 
to  use  program  PVRCK  and  evaluate  the  minimum  buckling  load  for 
each  value  of  axial  (m)  and  circumferential  (n)  Fourier  harmonic. 
Since  Donnell’s  equations  are  Inaccurate  for  n  _<  2,  no  attempt 
should  be  made  to  extrapolate  lower  than,  say,  n  =  4. 

The  selected  solution  procedure  was  to  determine,  for  each 

value  of  m  and  n,  the  minimum  buckling  mode  per  Eq.  (61), 

evaluate  b  per  Eq.  (73),  and  determine  X  from  Eq.  (75).  If  b  is 

s 

positive,  there  is  no  reduction  in  the  buckling  load  below  the 
classical  value.  This  procedure  was  followed  for  various  ranges 
of  Batdorf ' s  parameter  (.1  <  Z  <  1000)  and  imperfection 
geometries  (.0001  _<  E,  <  1.5). 

A  three-dimensional  perspective  plot  was  generated  for  each 
minimum  found.  Figure  2  illustrates  the  radical  behavior  of  b. 
For  different  values  of  Z,  the  behavior  appears  to  be  quite 
different.  For  Z  _<  2,  the  critical  wave  number  of  n  =  2  was 
determined.  Thus,  the  results  for  this  region  was  ignored. 

The  trends  of  Koiter’s  parameter  for  various  Batdorf's  para¬ 
meters  show  that  at  extremely  low  values  of  Z(<2),  b  is  either 
small  (>  -0.001)  or  positive.  Thus  the  shell  is  insensitive  to 
imperfections.  Some  moderate  fluctuations  occur  3  <  Z  <  10 . 

Note,  in  this  same  region  Hutchinson  predicted  the  greatest 
imperfection  sensitivity;  however,  this  was  never  confirmed  by 
experiment . 

For  low  values  of  Z  one  mode  appeared  to  dominate.  Above 
Z  >  20,  adjacent  modes  became  just  as  dominant  as  the  minimum 
mode.  For  higher  wave  numbers,  the  critical  mode  may  not  have 
the  largest  imperfection  sensitivity  factor. 

With  Increasing  Z,  the  predicted  critical  axial  mode  number 
also  increases,  verifying  Koiter’s  prediction  [11]: 
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Figure  2  Koiter's  Parameter  for  Axially  Loaded  Cylinder 


Figure  2  (continued) 


(e)  Z  =  20,  b  =  -  4.2 


(f)  Z  =  30,  b  =  -  10.98 


Figure  2  (continued) 


«  V 


U1 


(g)  Z  =  50,  b  =  -  20.09 


Figure  2  (continued) 


(i)  Z  =  100,  b  =  -  25.61 

Figure  2  (continued) 


(j)  Z  =  250,  b  =  -  4.82 


(JO  Z  =  500  ,  b  =  -  20.27 


U)  Z  =  750  ,  b  =  -  6.28 


Figure  2  (continued) 


(m)  Z  =  1000,  b  =  -  10.66 


Figure  2  (continued) 


✓z. 


m 


er 
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Since  these  results  are  based  upon  the  assumption  that 
boundary  conditions  have  little  influence,  it  is  anticipated  that 
the  predictions  are  more  conservative  than  those  obtained  with 
the  proper  boundary  conditions.  However,  the  magnitude  of 
Koiter's  parameter  which  corresponds  to  the  critical  mode 
(denoted  by  A  in  Pig.  3)  has  a  very  large  magnitude  and  is  quite 
erratic.  If  the  near  minimum  buckling  load  is  chosen  on  the 
basis  of  the  minimum  ordered  mode  (denoted  by  o  in  Pig.  3  which 
is  at  the  least  value  of  the  axial  wave  number),  then  the 
predictions  of  b  become  quite  well  behaved. 

This  same  observation  of  erratic  behavior  of  Koiter's 
parameter  was  reported  by  Yamaki  [19].  Two  different  analysis 
procedures  were  used  (an  asymptotic  and  a  full  nonlinear  method) 
and  the  results  obtained  were  similar  to  those  presented 
herein.  Yamaki 's  observation  was  that  the  lower  bound  of  the 
critical  loads  diminishes  significantly  with  the  increase  of  Z, 
as  well  as  R/t.  For  long  shells  with  Z  >  500,  the  imperfection 
sensitivity  factor  remains  constant.  Further,  the  postbuckling 
mode  is  always  symmetric  for  short  shells  Z  <  100;  while  for  long 
shells  with  Z  >  200  asymmetric  modes  are  predominant.  The  wave 
number  Immediately  after  buckling  is  generally  smaller  than  the 
critical  wave  number.  The  imperfection  sensitivity  parameter  for 
pure  torsion  and  hydrostatic  pressure  as  a  function  of  Z  was  well 
behaved,  while  for  compression  some  erratic  behavior  was 
observed . 
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Koiter's  Parameter 


Figure  3  Variation  of  Koiter's  Parameter  for  Axially  Loaded  Cylinder 


5.0  CONCLUSIONS  AND  OBSERVATIONS 

Pour  computer  program  procedures  have  been  assembled  that 
permit  a  convenient  and  economical  evaluation  of  the  imperfection 
sensitivity  of  cylindrical  shells.  When  all  of  the  imperfection 
data  is  known,  reasonable  agreement  to  experimental  results  is 
found  for  these  programs.  Because  of  the  basic  assumptions  behind 
each  program,  certain  caution  must  be  taken  In  attempting  to 
appiy  the  results  directly.  Similarly,  extreme  care  should  be 
taken  when  evaluating  any  experimental  data  because  not  all  of 
the  significant  data  may  be  reported. 

During  the  assembly  and  verification  of  the  computer 
programs,  a  variety  of  test  cases  were  run.  Prom  these 
solutions,  some  observations  on  the  factors  that  affect  behavior 
of  the  imperfection  sensitivity  of  shells  were  made.  Some  of 
these  observations  are  noted  as  follows: 

(1)  For  the  general  types  of  membrane  producing  loads 
(axial  compression,  torsion,  and  hydrostatic  pressure), 

compression  is  associated  with  the  greatest 
imperfection  sensitivity. 

(2)  When  torsion  is  combined  with  axial  compression,  even 
greater  Imperfection  sensitivity  is  observed. 

(3)  The  effects  of  boundary  conditions  may  mask 
imperfection  sensitivity,  particularly  for  axially 
stiffened  shells. 

(4)  Shells  with  outside  ring  stiffeners  tend  to  develop 
more  imperfection  sensitivity  than  those  with  internal 
ring  stiffeners. 

(5)  Two-way  stiffening  is  far  more  effective  than  single 
direction  stiffening. 

(6)  For  certain  cases,  a  multiplicity  of  eigenvalues  can 
exist  at  the  same  buckling  load.  For  those  states  of 
stress,  the  imperfection  sensitivity  corresponding  to 
these  adjacent  eigenvalues  can  be  considerably 
different . 
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(7)  No  logical  conclusions  can  be  made  for  one  shell 
geometry  which  are  applicable  to  all  types  of 
cylindrical  sizes  and  loading  conditions.  More 
extensive  parametric  studies  should  be  performed  in 
order  to  draw  general  conclusions. 

(B)  Par  more  detailed  experimental  evidence  is  required  to 
fully  demonstrate  a  particular  behavior  of  imperfection 
sensitivity.  The  effects  of  axisymmetric  imperfection 
have  been  fairly  adequately  documented.  The  effects  of 
asymmetric  imperfection  have  not. 

The  degradation  of  the  critical  axial  compression  load  not 
only  depends  upon  the  value  of  the  imperfection  sensitivity 
factor  but  also  on  the  magnitude  of  the  imperfection  and  its 
waveform.  Single  mode  behavior,  particularly  for  the  asymmetric 
mode,  suggests  that  only  small  imperfection  amplitudes  ought  to 
be  considered  when  employing  Koiter's  method.  Further  multi-mode 
participation  will  have  a  much  greater  affect  on  the  critical 
buckling  mode  than  single  imperfection  geometry. 

The  study  of  imperfection  sensitivity  of  shells  cannot  be 
adequately  resolved  by  simply  modifying  classical  linear  buckling 
theories  or  even  nonlinear  buckling  theories.  A  more  unified 
approach  of  handling  imperfection  geometry  is  required.  When  a 
complete  description  of  an  imperfection  geometry  is  established, 
adequate  correlation  to  experimental  test  results  can  be 
obtained.  There  is  no  current  accepted  method  that  can  be  used 
to  suggest,  in  advance,  which  critical  participating  modes  should 
be  used  in  the  analysis.  At  present,  this  can  be  accomplished 
only  through  trial  and  error. 
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